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[i] We summarize observations by the MESSENGER spacecraft of highly coherent waves 
at frequencies between 0.4 and 5 Elz in Mercury’s inner magnetosphere. This survey covers 
the time period from 24 March to 25 September 2011, or 2.1 Mercury years. These waves 
typically exhibit banded harmonic structure that drifts in frequency as the spacecraft 
traverses the magnetic equator. The waves are seen at all magnetic local times, but their 
observed rate of occurrence is much less on the dayside, at least in part the result of 
MESSENGER’S orbit. On the nightside, on average, wave power is maximum near the 
equator and decreases with increasing magnetic latitude, consistent with an equatorial 
source. When the spacecraft traverses the plasma sheet during its equatorial crossings, 
wave power is a factor of 2 larger than for equatorial crossings that do not cross the plasma 
sheet. The waves are highly transverse at large magnetic latitudes but are more 
compressional near the equator. However, at the equator the transverse component of these 
waves increases relative to the compressional component as the degree of polarization 
decreases. Also, there is a substantial minority of events that are transverse at all magnetic 
latitudes, including the equator. A few of these latter events could be interpreted as ion 
cyclotron waves. In general, the waves tend to be strongly linear and characterized by 
values of the ellipticity <0.3 and wave-normal angles peaked near 90°. Their maxima in 
wave power at the equator coupled with their narrow-band character suggests that these 
waves might be generated locally in loss cone plasma characterized by high values of the 
ratio (3 of plasma pressure to magnetic pressure. Presumably both electromagnetic ion 
cyclotron waves and electromagnetic ion Bernstein waves can be generated by ion loss 
cone distributions. If proton f3 decreases with increasing magnetic latitude along a field 
line, then electromagnetic ion Bernstein waves are predicted to transition from 
compressional to transverse, a pattern consistent with our observations. We hypothesize 
that these local instabilities can lead to enhanced ion precipitation and directly feed field- 
line resonances. 

Citation: Boardsen, S. A., J. A. Slavin, B. J. Anderson, H. Korth, D. Schriver, and S. C. Solomon (2012), Survey of coherent 
~1 Hz waves in Mercury’s inner magnetosphere from MESSENGER observations, J. Geophys. Res., 117, A00M05, 
doi: 10.1 029/20 1 2JA0 17822. 


1. Introduction 

[ 2 ] Highly coherent waves with frequencies near ~1 Hz 
were observed in Mercury’s magnetosphere during the first 
(Ml) [, Boardsen etal., 2009a] and second (M2) [, Boardsen etal. 
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2009b] near-equatorial flybys of the MErcury Surface, Space 
ENvironment, GEochemistry, and Ranging (MESSENGER) 
spacecraft. They were also briefly (~38 s burst) observed 
near an equatorial crossing during the first Mariner 10 flyby 
of Mercury in 1974 [ Russell , 1989]. During both Ml and 
M2, narrow-band, ultra-low-frequency (ULF) waves were 
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detected from just after closest approach nearly to the 
outbound magnetopause crossing. Bursts of waves were 
recorded almost continuously during Ml and M2 over 
approximately 9 min intervals on both encounters. During 
these wave observations MESSENGER covered a range in 
radial distance from 1.1 to 1.9 R M (where R M is Mercury’s 
radius), in magnetic local time (MLT) from 2 to 6 h, and in 
magnetic latitude from —23° to —15°. The wave frequency 
varied between the He + cyclotron frequency ,/ C iic+ and the 
proton cyclotron frequency f c H +. An outbound boundary 
layer was encountered on both flybys, and wave power was a 
factor of — 10 greater inside the boundary layer than outside 
the layer. No frequency drift was observed within the 
boundary layer, whereas a frequency drift was observed 
during Ml in the region inward of the boundary layer. The 
magnetosphere was quiet during Ml and highly disturbed 
during M2, and the overall wave power was 4 to 5 times 
larger during M2 than during Ml. Near the closest approach 
(CA) of each flyby trajectory, the power in the wave com- 
ponent oriented parallel to the local magnetic field tended 
to be larger than the power in the component perpendic- 
ular to the field, whereas away from CA the perpendicular 
power dominated. There was a large scatter in both the 
ellipticity and wave-normal angle, but the ellipticity values 
were skewed toward right-handedness, and the wave-normal 
angles tended to be greater than 45°. 

[ 3 ] There have been many explanations for these waves. 
Russell [1989] proposed that they are standing Alfven 
waves. However, due the substantial field-aligned fluctua- 
tions, Blomberg [1997] and Southwood [1997] argued 
against an explanation involving purely standing Alfven 
waves. Because these waves are primarily seen at frequen- 
cies between / cHe + and^ H +. and because a sizable heavy ion 
component has been documented in Mercury’s magneto- 
spheric plasma [Zurbuchen et al., 2008; Raines et al., 2011], 
the effects of heavy ions must be important in wave gener- 
ation. Zurbuchen et al. [2011] recently reported that heavy 
ions with values of mass per charge near that of Na , termed 
the Na + group, contribute substantially to the total plasma 
pressure in the near-Mercury nightside plasma sheet. 

[ 4 ] A number of studies have employed cold-plasma the- 
ory to explore the effects of ions on waves in the light-ion 
frequency range [Othmer et al., 1999; Glassmeier et al., 
2003; Kim and Lee, 2003; Klimashkin et al., 2006; Kim et 
al., 2008, 2011], Othmer et al. [1999] and Klimashkin et 
al. [2006] found that crossover frequencies, across which 
the polarization of a wave mode changes sign, are preferred 
frequencies for modified field-line resonances. However, 
Kim et al. [2008] showed from simulations that coupling 
occurs at ion-ion resonance frequencies. Kim et al. [2011] 
showed further that this ion-ion hybrid resonance coupling 
at Mercury is more efficient than the Alfven resonance for a 
plasma composed of H + and Na + , and that the field-aligned 
resonant frequency structure can be complex and discontin- 
uous and highly dependent on the density ratio of Na + to H + . 
An ion-ion hybrid frequency is located between each pair of 
neighboring ion cyclotron frequencies and just above the 
Bushsbaum-Bers frequency for that ion pair (a discussion of 
the difference between Buchsbaum-Bers resonances and 
ion-ion hybrid resonances is given by Lee et al. [2008]). If 
these waves are field-line resonances, their spatial distribution 


could be used to demarcate boundaries between open and 
closed field lines. 

[ 5 ] On the basis of MESSENGER flyby observations, 
Boardsen et al. [2009a, 2009b] suggested that because of 
their substantial compressional component these oscillations 
are not field-aligned resonances but rather quasi-trapped 
waves, similar to fast magnetosonic waves observed at fre- 
quencies near / cII+ near the equatorial plasmapause of Earth. 
Moreover, Boardsen et al. [2009a, 2009b] predicted that the 
waves would be confined in latitudinal extent. In Mercury’s 
inner magnetosphere the planetary loss cone half angle var- 
ies from 67° at an L-shell of 1.1 to 20° at an L-shell of 2.0 
(where the L-shell parameter defines an axisymmetric sur- 
face of those lines of magnetic force from the dipole com- 
ponent of Mercury’s internal field that intersect the magnetic 
equator at a distance L R M from the dipole center), leading to 
strongly non-Maxwellian velocity distributions that are 
expected to be highly unstable. Such anisotropic plasma 
distributions, with a temperature in the direction perpendic- 
ular to the local magnetic field greater than that in the par- 
allel direction, have been observed in the global hybrid 
simulations of the flybys [Trdvnlcek et al., 2010]. These 
instabilities lead in turn to the development of plasma waves 
that can scatter ions into the loss cone from which they 
precipitate onto Mercury’s surface [Schriver et al, 2011]. 

[6] That these ULF waves at Mercury are poorly under- 
stood is at least in part because of the limited spatial and 
temporal coverage of the flybys. The insertion of the MES- 
SENGER spacecraft into a near-polar orbit in March 2011 
has allowed mapping of these waves against magnetic lati- 
tude and magnetic local time. Here we characterize the 
spectral content, amplitudes, and polarization properties of 
these waves on the basis of orbital observations. 

[ 7 ] The variation of the MESSENGER orbit about Mercury 
as Mercury orbits the Sun is shown in Figure 1. The orbital 
plane of MESSENGER lies in the dawn-dusk plane when 
Mercury is near perihelion (0.3 1 AU) and aphelion (0.47 AU). 
Near perihelion, MESSENGER’S periapsis is on Mercury’s 
dawn side, and near aphelion MESSENGER’S apoapsis is on 
the dusk side. The orbital plane of MESSENGER lies near the 
noon-midnight plane of Mercury when the planet’s heliocen- 
tric distance is -0.37 AU. Near 0.37 AU MESSENGER 
samples either the near tail (with Mercury traveling outbound 
from the Sun) or the deep tail (with Mercury traveling toward 
the Sun) of Mercury’s magnetosphere. This information, 
combined with theory, aids in the physical interpretation of 
these waves, and thus in determining the role they play in 
Mercury’s magnetosphere. In this paper we present examples, 
statistics, and our interpretation of these waves from observa- 
tions during MESSENGER’S first six months in orbit. 

2. Data Analysis 

[8] We searched the full set of MESSENGER Magne- 
tometer [Anderson et ah, 2007] data from 3 March to 25 
September 2011, or 2.1 Mercury years, for plasma waves 
within Mercury’s magnetosphere at a radial distance of less 
than 2 R M from the planet’s center. The sampling rate of the 
Magnetometer data was 20 s -1 . Individual fast Fourier 
transform (FFT) spectra were generated for 20-s segments, 
over which the data had first been linearly detrended. Suc- 
cessive spectra were stepped 2.5 s in time from one to the 
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2011 091 (04/01) 0.37 AU 2011 119 (04/29) 0.47 AU 



201 1 146 (05/26) 0.37 AU 201 1 162 (06/1 1) 0.31 AU 



Figure 1. The MESSENGER orbit on days of the year in 
2011 when the spacecraft’s orbital plane lay approximately 
in the MSO (a) X-Z, (b) Y-Z, (c) X-Z, and (d) Y-Z planes as 
Mercury orbited the Sun. The orbital trajectory is in the 
counter-clockwise direction in Figures lb and lc and clock- 
wise in Figures la and Id. The radial distance of Mercury 
from the Sun is given in AU. Also shown are the intersec- 
tions of the (outer) bow shock and (inner) magnetopause 
from the model of Slavin et al. [2009]. 

next. For each time and frequency step, polarization para- 
meters (eigenvalues, eigenvectors, ellipticity, degree of 
polarization or coherency, and wave-normal angles) were 
computed from the average spectral power matrix. The 
method we used to compute these quantities is that described 
by Means [1972] and Arthur et al. [1976]. With this method, 
the real part of the spectral power matrix is diagonalized, and 
the wave-normal direction is given by the eigenvector of the 
minimum eigenvalue. With the rotation matrices that result 
from this diagonalization, the complex spectral matrix is 
then rotated into this diagonalized frame, and the submatrix 
given by the intermediate and maximum eigenvectors is 
used to compute the ellipticity and the degree of polarization 
(DOP). The average spectral power matrix for each fre- 
quency and time step is the average of the spectral power 
matrices over that particular frequency and time step and that 
of its two nearest neighbors in frequency and six nearest 
neighbors in time. For this averaging algorithm, we com- 
puted the coherency of a white-noise signal to be 0.41 with a 
standard deviation of 0.15. 

[ 9 ] The Magnetometer data are in Mercury solar orbital 
(MSO) coordinates, in which X is directed from the center of 
the planet toward the Sun, Z is normal to Mercury’s orbital 
plane and toward the north celestial pole, and Y is in the 
direction opposite to Mercury’s orbital motion. The alter- 
native Mercury solar magnetic (MSM) coordinate system is 


used for computing the magnetic equator, MLT, magnetic 
latitude, L-shell, and invariant latitude (cos' 1 ( l/s/L)). The 
magnetic dipole has been shown to lie nearly along the spin 
axis, and its location is offset from the planet center along 
the Z MSO axis by 0.19 R M [Anderson et al., 2008, 2010, 
2011]; i.e., the magnetic equatorial plane is parallel to the 
MSO X-Y plane but is offset by 0.19 R M in the +Z MSO 
direction. The transformation from MSO to MSM coordi- 
nates is thus simply 2f MSM = W MSO , Tmsm = Tmso> Zmsm = 
•Z’mso — 0.19 R m . 

[ 10 ] In addition to waves, MESSENGER frequently 
detected a constant tone near 0.67 Hz and its second har- 
monic, neither of which were observed during the flybys. 
This tone is the result of oscillations of the 3.6-m boom at 
the end which the Magnetometer sensor is mounted. For 
small-amplitude boom oscillations in a constant magnetic 
field, a transverse signal will be generated at the boom 
oscillation frequency. We examined 10 of these boom 
oscillation events, and all were associated with slew man- 
euvers about the spacecraft’s x-axis, which is aligned with 
the spacecraft solar panel arms and is perpendicular to the 
Magnetometer boom [Anderson et al, 2007]. Because the 
oscillations end abruptly at the end of a maneuver, the boom 
oscillations are believed to be fed by enhanced jitter from the 
spacecraft’s momentum wheels during the maneuver. The 
two hinge joints (necessary for deployment) of the boom are 
oriented in the spacecraft’s y-z plane, a location that renders 
the boom susceptible to boom oscillations during spacecraft 
slews about the x-axis. The choice of a 20-s interval for FFT 
analysis allows us to better separate boom oscillations from 
magnetospheric waves. Boom oscillations are not included 
further in our statistical study. 

3. Observations 

[ 11 ] An example of these waves as seen in dynamic 
spectra computed from the Magnetometer data is shown in 
Figure 2. On 29 March 2011, the spacecraft encountered 
narrow-band waves as it approached the magnetic equator 
from the north. Waves transverse to the local magnetic field 
were first seen at intermediate magnetic latitudes (starting at 
~35°N), and as the spacecraft approached the magnetic 
equator the wave amplitude increased and the waves shifted 
to compressional orientations (parallel to the local magnetic 
field) as the frequency shifted toward /c H+ . The fundamental 
mode and its second harmonic are clearly visible in Figure 2. 
During this and other wave observations, the fundamental 
frequency lay between f . He+ and f . H+ , and its frequency 
variation in time was typically smooth. The fundamental 
frequency often started near or below / cHe + at higher lati- 
tudes, rose to f cii+ as the spacecraft moved equatorward, 
and then tended to track f cii+ as the spacecraft moved over 
the southern hemisphere. Similar, highly coherent magnetic 
waves at frequencies near 1 Hz were nearly always observed 
(^80% of the equatorial crossings) within Mercury’s inner 
magnetosphere on the nightside, and generally these waves 
exhibited harmonic banded structure along with a frequency 
drift in the dynamic spectra. 

[ 12 ] A segment (50 s long) of the time series of the mag- 
netic field near the magnetic equator (^02:36 UTC) for this 
example of wave observations is shown in Figure 3. The 
data have been rotated into directions parallel and 
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HH:MM:SS 02:20:00 

02:25:00 

02:30:00 

02:35:00 

02:40:00 

02:45:00 

02:50:00 

L-SHELL 

13.66 

2.16 

1.28 

1.26 

1.52 

1.96 

237 

M-I AT 

74.79 

48.33 

23.55 

4.12 

-10.46 

-21.54 

-30.20 

MLT 

23.69 

1.22 

1.59 

1.78 

1.91 

2.01 

2.10 

R* 

0.94 

0.95 

1.07 

1.26 

1.47 

1.69 

1.92 

R 

1.12 

1.10 

1.16 

1.29 

1.45 

1.63 

1.83 


Figure 2. Dynamic spectra of the (top) parallel (compressional) and (middle) perpendicular (transverse) 
components and (bottom) degree of polarization (DOP) of the magnetic field during the observation of 
ULF waves at frequencies between 0.2 and 10 Hz following periapsis passage on 29 March 2012. The ver- 
tical black line indicates the magnetic equator. The white curves show the (top to bottom) H + , He ++ , and 
He + cyclotron frequencies. The time-axis label L-SHELL is the magnetic L-shell, M-LAT is magnetic lat- 
itude in degrees, MLT is magnetic local time in hours, R is the radial position of the spacecraft from planet 
center in R M , and R* is the radial distance from the offset dipole. Narrow-band harmonic waves are 
observed in all panels. The waves are highly compressional (see Figure 3) near the magnetic equator 
and transverse away from the magnetic equator. 


perpendicular to the average local magnetic field during this 
interval. Two-dimensional minimum variance analysis 
(MV A) has been performed on the two components in the 
perpendicular plane; Bj_i points in the minimum variance 
direction in that plane, and B , 2 points in the maximum 
variance direction in that plane. The compressional By 
component dominates, with a maximum peak-to-peak 
amplitude of ~10 nT in an ambient field of 84 nT, indicating 
that the wave magnetic field components are strongly linear 
and parallel to the local field. The compressional component 
is highly distorted from a sinusoidal waveform, consistent 
with the multiple harmonics observed in Figure 2. A mild 
magnetic depression marks the crossing of the plasma sheet 
near the equatorial crossing. 

[ 13 ] For a substantial minority (~25%) of the wave 
events, the compressional component was not dominant near 
the magnetic equator but instead the transverse component 
was larger. An example of such an event on 1 August 2011 
is shown in Figure 4. During this event the transverse com- 
ponent dominated throughout, from 09:09 to 09:25 UTC. A 
time series of the largest-amplitude waves, about ^6 nT 
peak-to-peak oscillation in a 90 nT ambient field, is shown 
in Figure 5. The near sinusoidal signal observed in all three 
components makes this an ideal event for MV A. A hodo- 
gram of the intermediate ( B 2 ) and maximum (/i 3 ) variance 
components of this time series is given in the inset in 
Figure 5. The ratio of intennediate to minimum eigenvalue is 
12, which yields a well-detennined wave-nonnal angle 


(WNA) of 35°. The waves at the magnetic equator are left 
handed with respect to the ambient magnetic field and have 
an ellipticity of —0.35, which is suggestive of the ion 
cyclotron mode. However, MVA applied to a few wave 
cycles at a time away from the equator yielded mixed results, 
with both left- and right-handed polarization. Near the 
equator, the orientation of the major axis is perpendicular to 
our estimate of the gradient in the magnetic field magnitude. 
However, this example is an exception in that even events 
that are dominantly transverse at the equator tend to be 
highly linear. For the data set analyzed for this paper, only 
two other orbits were identified with strong transverse 
oscillations at the equator; for both such events the major 
axis of the polarization ellipse was also perpendicular to our 
estimate of the gradient of the magnetic field magnitude. 

[ 14 ] We did not clearly detect waves in the tail at distances 
R > 2 R m , in settings for which the orbit was similar to that 
of Figure lc. A few examples of two to four cycles of low- 
amplitude oscillations at a frequency near 1 Hz were 
observed in the tail at such distances. A statistical analysis of 
the distribution of observed waves, presented next, is fol- 
lowed later in the paper by a brief discussion of possible 
reasons why these waves are not seen deeper in the tail. 

4. Statistics 

[ 15 ] The overall coverage of the wave observations is 
shown in Figure 6. Average integrated power spectral 
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HH:MM:SS 02:36:19 

02:36:32 

02:36:44 

02:36:57 

02:37:09 

Xmso(Rm) 

-1.17 

-1.17 

-1.18 

-1.19 

-1.19 

Ymso(Rm) 

-0.60 

-0.61 

-0.61 

-0.62 

-0.62 

Zmso(Rm) 

0.19 

0.17 

0.16 

0.14 

0.13 

R(Rm) 

1.33 

1.33 

1.34 

1.35 

1.35 


Figure 3. Strongly compressional waves near the magnetic equator for a portion of the data set depicted in 
Figure 2. The first and second panels are the perpendicular components of the field (see text), the third panel 
is the parallel component, and the fourth panel is the field magnitude. Peak-to-peak amplitudes are ^ 1 0 nT. 


density between 0.8 and 5.0 Hz binned by magnetic latitude 
and MLT is shown in Figure 6a. Because of MESSENGER’S 
orbit and the proximity of Mercury’s dayside magnetopause 
to the planet, there was no Magnetometer coverage near the 
dayside magnetic equator. The plot shows clear enhancements 
in intensity in two regions. One is more or less centered about 
the magnetic equator, the result of these narrow-band ULF 
waves, and is concentrated on the nightside (18 to 6 h MLT). 
The second is in the cusp region (9 to 14 h MLT at magnetic 
latitude >40°). From visual inspection of the dynamic spec- 
trum in the cusp region, the magnetic fluctuations in the cusp 
region tend to be broadband. 

[i6] We visually examined each dynamic spectrogram for 
each orbit in the data set, and we marked regions in time and 
frequency at which narrow-band waves were observable, 
taking care to exclude boom oscillations. Within each 
marked region the following criteria were used to identify 
wave events: (1) DOP >0.90 (about three standard devia- 
tions above the white-noise level); (2) the power spectral 
density >0.2 nT 2 /Hz; and (3) the power spectral density was 
more than two standard deviations above the estimated 
background of the spectrum. Of the frequencies at which 
these criteria were satisfied, the frequency at maximum 
power was retained for each time step. With these criteria 
about 21,000 wave events were identified. The distribution 
of those events in MLT and magnetic latitude is shown in 
Figure 6b. The waves were detected at all MLTs, but they 
were concentrated on the nightside. For the wave events near 
midnight MLT the orbital configuration was similar to that 
shown in Figure la, and the waves occurred mainly in 
northern hemisphere regions where the magnetic field lines 


were most likely closed. From an MLT of 20 h to an MLT of 
14 h in Figure 6b the distribution of events shows a tendency 
to split into two branches, one near the magnetic equator and 
another that rises to higher latitudes from dusk to noon. This 
apparent splitting was seen in both Mercury years. A similar 
splitting is seen between MLT noon and dusk. 

[n] The location of plasma sheet crossings during this 
period for R M < 2 is indicated in Figure 6c. As expected 
these crossings are located on the nightside, and the highest 
concentration is seen around midnight. Plasma sheet inter- 
vals are identified as regions in which the magnetic field 
strength is at least 10 nT less than the adjacent tail-lobe field 
strength. The plasma sheet was encountered only on about 
half of the nightside crossings of the magnetic equator. A 
similar distribution of plasma sheet locations was given by 
Korth et al. [2011], 

[is] Histograms of the number of wave events as a func- 
tion of MLT, magnetic latitude, and radial distance are given 
in Figures 7a-7c. For normalization, the number of time 
samples in the data set is given as a function of MLT and 
magnetic latitude in Figures 7d-7f. In Figures 7g-7i the 
number of events (Figures 7a-7c) is shown normalized by 
respective coverage (Figures 7d-7f). The distribution of 
wave events with MLT drops off strongly on the dayside, 
but this pattern is at least in part the result of the lack of 
coverage near the dayside magnetic equator as well as 
broadband turbulence in the cusp region masking these 
waves if they occur. The distribution of events peaks near 
the magnetic equator (magnetic latitude 10 ± 18°). The 
distribution of events peaks at a radial distance of (1.35 ± 
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Figure 4. An example of a ULF wave for which the transverse component was dominant at all latitudes, 
even during equatorial passage. In addition to the magnetospheric waves, oscillations of the Magnetometer 
boom induced by spacecraft attitude adjustments are visible at a frequency of ~0.7 Hz at 9:08 UTC. 
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Figure 5. Strongly transverse waves near the magnetic equator for a portion of the data set depicted in 
Figure 4. Maximum peak-to-peak amplitudes are ~6 nT. Inset: Hodogram of the intermediate (B 2 ) and 
maximum ( // ) variance components of the time series. The major axis of the polarization ellipse is per- 
pendicular to the estimated gradient in field magnitude. The wave frequency /is ~1.1 Hz, the ratio // m 
is ~0.8, the WNA is ~35°, and the waves are left handed with an ellipticity of —0.35. 
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Figure 6. Coverage in magnetic latitude and magnetic local time from 24 March-25 September 2011, for 
periods when the spacecraft was inside Mercury’s magnetosphere and R <2 R M . (a) Integrated spectral 
power over frequencies between 0.8 and 5.0 Hz for all data during this period (both with and without 
waves), (b) Observed wave events versus magnetic latitude and magnetic local time. Note the bifurcation 
in latitude of the distribution between dawn and noon and again from 20 h MLT toward dusk, (c) Plasma 
sheet crossings versus magnetic latitude and magnetic local time. 


0.13) R m - The frequency of events drops toward zero as the 
radial distance increases toward 2 R M . 

[ 19 ] Histograms of several characteristics of the wave 
events are given in Figure 8. A histogram of wave frequency 
(Figure 8a) shows three distinct peaks: one at a frequency 
just below f c ne++, a dominant peak at a frequency just below 
f c n+, and a peak near the second harmonic of that frequency 
at ~2/ C | |4 . A histogram of the ratio of power in the parallel 
component to the total power (Figure 8b) indicates that the 
number of compressional events approximately equals the 
number of transverse events. If the DOP selection criterion is 
raised above 0.90, the distribution of events becomes 
skewed toward the compressional side (1.0), and if the DOP 
selection criterion is reduced below 0.90 the distribution of 
events becomes skewed toward the transverse side (0.0). A 
histogram of wave ellipticity (Figure 8c) indicates that the 
waves tend to be linear, with ellipticity peaked near zero 
(0.10 ± 0.37) and with some bias toward right-handed 
polarization. A histogram of wave-normal angles (Figure 8d) 
shows a peak near 90° . There is also a small peak near ~ 10° . 
For transverse events that are strongly linear, the error in the 
wave-normal angle is large and contributes to the broad 
distribution observed at angles between 0° and 75°. 

[ 20 ] One of the most dominant features of the distribution 
of wave events is that the waves are strongly compressional 
near the magnetic equator and strongly transverse far from 
the equator. This pattern can be seen in a plot of the ratio of 
parallel-component power to total power for nightside 
events (Figure 9a). Near 0° latitude, most values of this 
power ratio exceed 0.7; as the DOP selection criterion is 
raised this tendency is strengthened, whereas as the DOP 


selection criterion is lowered this tendency is weakened. The 
total power in the wave events also maximizes around the 
equator and drops off with increasing magnetic latitude, as 
shown in Figure 9b. 

[ 21 ] The average wave frequency and magnetic field 
strength for the 969 events inside the plasma sheet is 0.90 ± 
0.36 Hz and 59 ± 13 nT, respectively. In contrast, the aver- 
age wave frequency and magnetic field strength of the 
8625 nightside wave events within 10° of the magnetic 
equator but on orbits that did not encounter the plasma sheet 
is 1 .27 ± 0.6 Hz and 86 ± 20 nT, respectively. The average 
power spectral density in the plasma sheet is about a factor of 
2 larger than that outside at similar latitudes. This finding is 
qualitatively consistent with that of the flybys. The average 
ratio of the compressional to total wave power is 0.6 near the 
equator for both the flybys and orbital crossings. 

[ 22 ] A plot of ellipticity versus magnetic latitude for wave 
events that are strongly transverse (ratio of power in the 
perpendicular component to total power >0.7) is shown in 
Figure 9c. Near the magnetic equator the absolute value of 
the ellipticity is ~0.5, whereas at higher magnetic latitudes 
the absolute value of the ellipticity decreases toward 0 as the 
waves become more strongly linear. Shown in Figure 9d is a 
histogram of the angle (() between the maximum eigen- 
vector direction and the estimated direction of the magnetic 
field gradient for those events for which the angle between 
the maximum eigenvector direction and the local magnetic 
field direction is greater than 70° and the magnetic latitude is 
greater than 20°. The median of ( is 70°, a value suggesting 
that far from the magnetic equator the maximum wave 
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Distribution of Wave Events 



Coverage Inside Magnetosphere R < 2 R M 



Figure 7. Histograms of wave events as a function of (a) magnetic local time, (b) magnetic latitude, and (c) 
radial distance; coverage in (d) MLT, (e) magnetic latitude, and (f) radial distance; (g-i) wave events from 
Figures 7a-7c nonnalized by their respective coverage from Figures 7d-7f. The total number of time samples 
with wave events is 21,482 for Figures 7a-7c, and the total number of samples is 381,981 for Figures 7d-7f. 
The percent occurrences in Figures 7g-7i correspond to bin sizes of 0.5 h, 2.5° latitude, and 0.02 R M radial 
distance, respectively. In Figures 7b and 7c median, mean, and standard deviation (Stdv) are also given. 


magnetic field fluctuations are nearly perpendicular to the 
background magnetic field gradient. 

5. Discussion 

[ 23 ] The plasma and magnetic field properties of Mer- 
cury’s inner magnetosphere provide important context to the 
wave observations reported in this paper. Mercury occupies 
a large fraction of the volume of this region; scaled to 
Earth’s magnetosphere the surface of Mercury would be at 
radial distance of 8 R E [Slavin et al., 2008], where R E is 
Earth’s radius. This large volume fraction filled by the planet 
means that on the dayside there is no minimum in magnetic 
field magnitude at the magnetic equator [Alexeev et al., 
2008] and that on the nightside the true magnetic equator 
(the point along a field line at which the background mag- 
netic field strength is a minimum) can have large fluctua- 
tions about the mean position in response to changing solar 
wind conditions. Because compressional waves tend to be 
confined near magnetic field minima [see Boards en et al., 
1992, Figures 4-6], they should not be observed at low 
latitudes on the dayside if generated locally. Also the tran- 
sition from compressional to transverse (Figure 9a) waves 
with increasing magnetic latitude can be understood at least 
qualitatively in terms of the confinement of compressional 


waves about the magnetic field minimum. Because the 
plasma at Mercury has multiple ion species, as transverse- 
dominant plasma waves (Figure 9c) move away from the 
equator into regions of increasing field strength the ratio of 
wave frequency to f cli+ decreases and waves approach a 
plasma crossover frequency (e.g., between / cH4 and f C He++, 
or between /[. H e++ and f c Na +). As waves approach such a 
crossover frequency they become more linear [Garnett et al., 
1965; Rauch and Roux, 1982; Stix, 1992]. This explanation 
at least qualitatively allows one to understand the trend in 
the absolute value of the ellipticity (Figure 9c), which 
decreases toward zero (linear waves) with distance from the 
magnetic equator. Electromagnetic ion cyclotron waves 
(ICWs), which have ellipticity <0, should be guided along a 
magnetic field line [Rauch and Roux, 1982], but for the 
results in this paper there is a skew toward right-handedness, 
even for transverse-dominant waves (Figure 9b). Unless 
there are field-aligned mass density irregularities present, 
these waves will not follow the field line. For a source that is 
distributed in space for non-ducted waves, wave spectra 
would be broadband away from the source, in contrast to the 
narrow-band spectra that are observed. We note, however, 
that the observation of right-handed transverse waves could 
be due to the limitations of the standard polarization analysis 
used in this study. 
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Figure 8. Histograms of wave events as a function of (a) frequency, normalized by the proton cyclotron 
frequency, (b) ratio of the parallel power to total power, (c) ellipticity (— 1 is left-handed circular, 0 is lin- 
ear, and 1 is right-handed circular polarized with respect to the background field), and (d) wave-normal 
angle. A total of 21,482 samples was used for each histogram. 


[ 24 ] Polarization analysis works best when one dominant 
wave mode is present. When two or more modes are com- 
peting for dominance, their superposition renders difficult 
the interpretation of wave-normal angle and ellipticity 
derived from such analysis [. Anderson etal., 1996; Denton et 
al., 1996]. Denton et al. [1996] investigated the case when 
two ICW modes are present with identical frequencies, 
amplitudes, ellipticities, and wave-nonnal angles. For this 
case the only free parameters are the angle 6a between the 
minor axes of the two polarization ellipses and the relative 
phase 6p between the two modes. Denton et al. [1996] 
showed that, as these two parameters are varied, wide ran- 
ges in wave-normal angle and ellipticity can result from the 
polarization analysis of the superposed waves, even for 
solutions with positive ellipticity. 

[ 25 ] If the transverse waves reported here are field-aligned 
resonances (i.e., standing waves constructed from ICWs), 
they can be regarded as a superposition of two ICWs prop- 
agating in opposite directions, similar to the model of 
Denton et al. [1996]. For ICWs with a wave-nonnal angle of 
0° (i.e., purely left-circularly polarized), their superposition 
will always yield left-circularly polarized waves. For non- 
zero wave-normal angles, however, part of (6a, 6<p) phase 
space will yield right-handed solutions from polarization 
analysis. We repeated the two-wave superposition analysis 
of Denton et al. [1996] using as input a wave-normal angle 
of 35° and an ellipticity —0.35 (from the left-handed waves 
shown in the inset of Figure 5). In that analysis 84% (16%) 
of (6a, 6ip) phase space yielded negative (positive) ellipticity 
values with a normalized median power of 0.68 (0.32). We 
compared this outcome with the data points in Figure 9c for 


magnetic latitudes between 0° and 15°, where there is a clear 
bifurcation in the ellipticity values. Of these points, 362 have 
negative ellipticity with a median power of 0.38 nT 2 /Hz, and 
1099 have positive ellipticity with a median power of 0.25 
nT 2 /Hz, i.e., 25% (75%) of the data points have negative 
(positive) ellipticities with a normalized median power of 
0.60 (0.40). 

[ 26 ] The median power indicated by polarization analysis 
of both the observations and the results from the two-wave 
superposition model is larger for left-handed than right- 
handed ellipticities. This difference can be understood as the 
superposition of two left-handed wave modes, each of which 
can be regarded as a superposition of a larger-amplitude left- 
handed circular trace and a smaller-amplitude right-handed 
circular trace. So, in order to add two such left-handed wave 
modes to obtain a right- (left-) handed solution from the 
polarization analysis of the sum of two left-handed waves, 
the two larger (smaller) left- (right-) handed traces must 
nearly cancel out to yield a smaller- (larger-) amplitude 
solution. A total of 84% of the phase space was left handed 
in the analysis of superposed model waves, compared with 
only 25% of the data points. This difference might be 
attributable to non-uniform sampling of (6a, 6ip) phase 
space, with both 6a and 6 ip constrained by the reflection of 
these waves off Mercury or because of approach to an ion- 
ion resonance along the field line. Although more work is 
needed on this topic, the transverse waves reported here may 
be primarily left handed rather than a mixture of right- and 
left-handed waves as indicated by the results of standard 
polarization analysis. 
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Figure 9. (a) Ratio of parallel power to total power versus magnetic latitude for nightside wave 
events. The waves tend to be compressional near the equator and strongly transverse away from the equator, 
(b) Average power spectral density of nightside wave events versus magnetic latitude, distinguished on the 
basis of position relative to the plasma sheet, (c) Wave ellipticity versus magnetic latitude for wave events 
for which the ratio of transverse power to total power is greater than 0.7. (d) Histogram of the angle (£) 
between the direction of the maximum eigenvector and the estimated direction of the gradient in the local 
magnetic field, for those events for which the angle between the maximum eigenvector direction and the 
direction of the local magnetic field is greater than 70° and the magnetic latitude is greater than 20°. 


[ 27 ] The observation that, far from the magnetic equator, 
the direction of the maximum eigenvector for wave magnetic 
field variation is nearly perpendicular to the background 
magnetic field gradient (Figure 9d) is strongly suggestive that 
there is a substantial component of field-aligned resonances 
to these wave observations. At frequencies below the lowest 
cyclotron frequency in the Alfven resonator frequency range, 
simulations predict that at the point where compressional 
waves couple with transverse waves the transverse magnetic 
field component of the wave will be strongly perpendicular to 
the background magnetic field gradient [Zhu and Kivelson, 
1989]. This behavior is also predicted for mode coupling at 
the ion-ion hybrid resonance [Kim et al., 2008]. 

[ 28 ] In a simulation with a constant background magnetic 
field, Kim et al. [2008, Figure 4] showed that compressional 
waves are mode converted to transverse waves for which the 
dominant transverse magnetic field component is perpen- 
dicular to the mass density gradient. The coupling occurs at 
the H + -Na + hybrid resonance at a frequency of 0.2 f c h+. 
Moreover, Kim et al. [2011, Figure lb] detennined the H + - 
Na hybrid resonance and its first two harmonics versus 
fractional composition of Na + in a H + -Na + plasma. From the 
theory of Kim et al. [2011] we may estimate the Na + frac- 
tional composition for the event shown in Figures 4-5. At 
the magnetic equator the observed flf C H+ ratio is ^0.8; at that 
value of the ratio in Kim et al. [2011, Figure lb], their theory 
predicts that the fractional composition of Na + is 0.8, 0.9, 


and 0.95 for the first three harmonics, respectively. A few 
estimates of the Na + ratio have been determined in Mer- 
cury’s near magnetosphere from measurements made by 
MESSENGER’S Fast Imaging Plasma Spectrometer (FIPS) 
during the Mercury flybys [Raines et al., 2011] and from 
orbit [Zurbuchen et al., 2011]. These estimates, summarized 
in Table 1, span the range 0.08 to 0.52, values lower than the 
values of 0.8 to 0.95 predicted by the Kim et al. [2011] 
simulations. The median value of the ratio f/f c n+ for all 
wave events analyzed here is 0.72 (Figure 8a), a value that 
suggests a larger fractional composition of Na + than what 
has been reported from analysis of the plasma ion data to 
date. However, an additional ion species added to the 
plasma, e.g., He ++ , would modify the ion-ion resonances, 
e.g., the H-He ++ and He ++ -Na + hybrid resonances would 


Table 1 . FIPS Observations in Mercury’s Plasma Sheet a 


Flyby 

« H +(cm 3 ) 

Th + (eV) 

% Na + 

Reference 

Ml 

1 tolO 

170 eV 

15 

Raines et al. [20 1 1 ] 

M2 

4 to 5 

690 eV 

8 

Raines et al. [20 1 1 ] 

In orbit 

1.6 and 0.28 

310 and 670 

1 1 and 52 

Zurbuchen et al. [2011] 


a The quantity n H+ is the proton density, 7^+ is the proton temperature, 
and % Na + is the fractional contribution to the plasma pressure by Na + - 
group ions and is derived under the assumptions that there is pressure 
balance between the plasma sheet and adjacent lobes and that Th+ ~ T Na+ 
[Raines et at., 201 1], 
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Figure 10. (a) Wavelength A, (b) ellipticity, and (c) ratio of the wave compressional field (<55ii) to the 
total wavefield ( SB ) versus frequency for a wave-normal angle of 35° in a background field of 90 nT. 
The plasma density is 1 cm -3 and the fractional ion composition is 0.86 H + , 0.05 He ++ , and 0.09 Na + . 
The crossover frequencies are indicated along with^ He + and^He-w- Blue denotes left-handed polarization, 
and red denotes right-handed polarization. 


replace the H + -Na hybrid resonance. For wave frequencies 
observed between f . He ++ and f cli + one would predict the 
fractional composition of He ++ as opposed to that of Na + . 

[ 29 ] If these transverse waves are field-line resonances 
composed of ICWs, how do the wavelengths compare with 
the length of the field line? From cold-plasma theory we 
have estimated the wavelengths, ellipticities, and compres- 
sional field in Figure 10. The wave-normal angle and 
background field are from Figure 5, but plasma density and 
ion composition are not known, so we adopted 1 cm~ 3 for 
the plasma density and assumed that the plasma is composed 
of H + , He ++ , and Na + . The ratio of He ++ to H + is taken to be 
5%, a typical solar wind value, and the ratio of Na to total 
ions is assumed to be 9% (see Table 1). The observed wave 
frequency in Figure 5 is 1.1 Hz, which gives a wavelength 
for the left-handed mode in Figure 10 of 0.34 R M and an 
eccentricity of —0.75. The approximate condition for a field- 
line resonance is 2/ DIP /A|| = n , where n is an integer, / DIP = 
1.8 R m is the length of the dipole field line, and Ay is the 
parallel wavelength. For a wavelength of 0.34 R M the field- 
line resonance harmonic number n is ~8, which is large. 
Note that the observed ellipticity of —0.35 from Figure 4 is 
smaller than —0.75 from Figure 10. A number of parameters 
can be changed to lower the magnitude of the model ellip- 
ticity, such as increasing the wave-normal angle, He ++ con- 
centration, or Na concentration. For demonstration purposes 
we do the last. When the Na + concentration is raised to 55% 
and the ratio of He ++ to H + is held fixed, the ellipticity 
becomes —0.37, the wavelength increases to 1.77 R M , and 
the harmonic number n of the resonance is ~2. Therefore it is 
plausible that these transverse waves could be field-line 
resonances. 

[ 30 ] The above arguments are based on cold-plasma the- 
ory, which is valid only for low values of the ratio (3 of 


plasma pressure to magnetic pressure. Because Mercury has 
no ionosphere and a low rotation rate, it is not thought to have 
a cold-plasma component inasmuch as there is no co- 
rotational electric field required to trap cold ions and no 
ionosphere that could serve as a source of low-energy ions. 
The ion plasma parameters at Mercury are important in 
understanding the plasma wave phenomena. Because of the 
limited field of view of the sensor, however, plasma para- 
meters have been estimated from FIPS measurements only 
for a limited number of intervals when MESSENGER was in 
Mercury’s plasma sheet during the Mercury flybys [Raines 
et al., 2011] or in orbit [Zurbuchen et al., 2011]; see Table 1. 
A determination of the relative composition H + and Na + is 
derived from the assumption of pressure balance between the 
magnetic pressure and particle pressure in the plasma sheet 
and lobes. From these parameters as well as typical values of 
the magnetic field magnitude B observed in the plasma 
sheet, we show in Table 2 estimates of the proton beta /3 H + 
derived from the ratio of estimated proton thermal speed Ch+ 
to the Alfven speed V A and the relation /3 H+ = (Ch+YFa) 2 . 

[ 31 ] A number of factors could explain the drop off in the 
occurrence of wave events with radial distance (Figure 7c). 


Table 2. Estimates of /3 H + for Several Combinations of Plasma 
Density n e and Fractional Ion Composition From Table 1 


B 

(nT) 

n e 

(cm" 3 ) 

% H + 

% Na + 

Va 

(km/s) 

C H+ (km/s) 

C h+ /T a 

/3h+ 

60 

5 

100 

0 

590 

170-250 

0.3-0.4 

0.09-0.16 

60 

5 

90 

10 

330 

170-250 

0.5-0. 8 

0.25-0.64 

60 

5 

50 

50 

170 

170-250 

1.0- 1.4 

1. 0-2.0 

60 

1 

100 

0 

1320 

170-250 

0.13-0.19 

0.02-0.04 

60 

1 

90 

10 

740 

170-250 

0.22-0.34 

0.05-0.12 

60 

1 

50 

50 

380 

170-250 

0.45-0.67 

0.2-0.45 
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HH:MM:SS 17:01:60 

17:06:60 

17:11:60 

17:16:60 

L-SHELL 

1.48 

1.27 

1.43 

1.77 

M-LAT 

31.99 

11.65 

-3.98 

-15.95 

MLT 

21.00 

21.22 

21.36 

21.47 

R* 

1.06 

~1 .22 

1.42 

1.63 

R 

1.17 

1.27 

1.42 

1.59 


Figure 11. A spectrogram (in the same format as Figures 2 and 4) that displays two intensity enhance- 
ments, mainly compressional and dominantly in the first two harmonics, which straddle the magnetic 
equator (black vertical line). The ratio of the fundamental frequency tof c n is ~0.7. 


If these waves are field-line resonances they will occur on 
closed field lines and the number of closed field lines will 
decrease with increasing radial distance down the tail. Near 
midnight these waves are not observed above an invariant 
latitude of 40°, which could indicate a boundary between 
open and closed field lines on the nightside. Another 
explanation is that a substantial number of protons are in 
bounce resonance with the wave period. In a dipolar field, 
500 eV protons will be in bounce resonance with the equa- 
torial proton cyclotron period at an L-shell of 3. If the proton 
temperature is 500 eV, there will be a substantial number of 
protons near 2 keV energy that will be in bounce resonance 
with the equatorial cyclotron period at an L-shell of 2. If 
these waves are generated locally by a loss cone instability, 
then from an L of 1 . 1 to 2 the planetary loss cone angle will 
decrease from 67° to 20°, providing a free energy source for 
the instability that decreases with increasing L-shell [e.g., 
Schriver et al., 2011]. 

[ 32 ] The mean radial location of these waves of 1.35 R M 
(Figure 7c) corresponds to a distance of 10.7 R E in Earth’s 
nightside magnetosphere. Narrow-band waves near the pro- 
ton cyclotron frequency have recently been reported in 
Earth’s plasma sheet boundary layer \Engebretson et al, 
2010] at distances of 8 to 18 R E in the tail, and both beam 
and ring (shell) proton velocity distributions were observed 
during these events. The frequency-time structure of these 
spectra is similar to that for the waves observed at Mercury, 
except that the waves at Earth are transverse dominant. From 
an instability analysis of these waves, Denton et al. [2010] 
concluded that they are generated by proton ring (shell) 
distributions for large /?h+- For /? H + ~ 0.4 the waves are 
electromagnetic transverse dominant, and as /3 H + increases 
they remain electromagnetic but become compressional 


when /? H + ~ 2 [Denton et al., 2010]. Gary et al. [2010] 
pointed out that this wave mode is the electromagnetic ion 
Bernstein mode. Whereas Denton et al. [2010] adopted an 
anisotropic ring distribution, Gary et al. [2010] modeled the 
same instability with an isotropic ring distribution for /3 H + ~ 
0.5 and obtained nearly identical growth rates to those of 
Denton et al. [2010], suggesting that anisotropy is not 
important for this instability. Gary et al. [2010] did not 
explore the instabilities at j3 a+ ~ 2, at which the waves 
became compressionally dominant [Denton et al., 2010]. In 
both studies ff c n+ at maximum growth rate was very close to 
1 and the wave-normal angle was near 90°; Gary et al. 
[2010] computed a wave-normal angle of 86° at maximum 
growth rate. Denton et al. [2010] found that these waves are 
highly linear and cited a typical ellipticity of 0.02. In our 
statistical study the median values from Figure 8 are 0.72 for 
flf c h+, 0.08 for ellipticity, and 7 1 ° for the wave-normal angle. 
Given the limitations of polarization analysis these values 
are consistent with the electromagnetic ion Bernstein mode. 
Gary et al. [2010] showed further that the fundamental mode 
does not always dominate. In our study the fundamental 
mode is usually but not always the dominant harmonic; in 
Figure 11, for instance, the power in the 2nd harmonic is 
larger than that of the fundamental mode just prior to the 
magnetic equator crossing. An approximate estimate of the 
wavelength A of the fundamental mode can be made from 
the relation [Gary et al., 2010, Figure 4] A ~ L|/(3/ cH+ ), 
where V\ is the ring velocity. For a ring velocity equal to a 
proton thermal speed of 210 km/s and a 60 nT field, the 
wavelength is ^80 km. If these waves correspond to the 
electromagnetic ion Bernstein mode, their wavelengths are 
much smaller than the dimension of Mercury’s magnetosphere. 
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[ 33 ] For the waves represented in Figure 9a the average 
background magnetic field energy density increases with 
increasing magnetic latitude. If the proton energy density is 
constant on average over this range of magnetic latitude, 
then /3 h + will decrease with increasing magnetic latitude. If 
these waves are mainly composed of electromagnetic ion 
Bernstein modes, then as /? H+ decreases they will transition 
from compressional to transverse dominant [Denton et al., 
2010]. Therefore, the transition from compressional to 
transverse dominance as magnetic latitude increases 
(Figure 9a) can be qualitatively understood if these waves 
are electromagnetic ion Bernstein modes in a high-/3 H + 
plasma for which /3 H + decreases with increasing magnetic 
latitude. Further experimental and theoretical analysis is 
warranted to explore this possible interpretation. 

[ 34 ] Shell distributions observed deep in Earth’s plasma 
sheet boundary layer (as opposed to the shell distributions 
observed near and inside geosynchronous orbit at Earth 
[Mcllwain, 1972] for which a co-rotational electric field is 
required, unlike the situation at Mercury) are believed to be 
generated by counterstreaming protons that are scattered in 
pitch angle [Engebretson et al, 2010]. Observations of such 
distributions at Mercury have not been reported, and 
whether a similar mechanism occurs in Mercury’s plasma 
sheet boundary layer that would generate proton shell dis- 
tributions has not been explored. However, the strong loss 
cone distributions predicted to occur in Mercury’s inner 
magnetosphere could provide the positive perpendicular 
slope necessary for the generation of waves for this /3 H +- 
dependent instability. Figure 1 1 shows a spectrogram for an 
event characterized by two intensity enhancements that are 
mainly compressional and are observed to straddle the 
magnetic equator. These events could be interpreted as 
waves generated by local instabilities during two crossings 
of the plasma sheet boundary layer, although the location of 
the plasma sheet and its boundary layers are not clear from 
only the Magnetometer data for this event. 

[ 35 ] The observation that the wave power is maximum at 
the equator and is a factor 2 larger in the plasma sheet sug- 
gests that the waves reported in this study are due to a local 
plasma instability. The increase in power at the equator 
could be attributed to the focusing of compressional waves 
from the magnetopause, but we believe that this focusing 
would produce a broadband spectrum rather than the nar- 
row-band spectra that are observed. If field-line resonances 
exist, they could be driven by local instabilities, in contrast 
to compressional waves that are generated remotely at the 
magnetopause. Linear stability analysis using a loss cone 
distribution in high-/? H + plasma is required to determine if 
loss cone distributions can generate both transverse and 
compressional electromagnetic waves with changing /? H +. 

6. Conclusions 

[36] This paper has presented a survey of highly coherent 
ULF waves at frequencies between 0.4 and 5 Hz in Mercury’s 
inner magnetosphere (R< 2 f? M )- The survey covers the time 
period from 24 March to 25 September 20 1 1 , or 2. 1 Mercury 
years. These waves usually exhibit banded harmonic struc- 
ture that drifts in frequency as the spacecraft traverses the 
magnetic equator. These waves are observed at all MLTs, but 
their occurrence strongly drops off on the dayside, which is 


partly due to orbital bias. On the nightside, on average, the 
wave power is maximum near the equator and decreases with 
increasing magnetic latitude, suggesting an equatorial source. 
When the spacecraft traverses the plasma sheet during its 
equatorial crossings, the wave power is a factor of 2 larger 
than equatorial crossings that do not cross the plasma sheet. 
The waves are highly transverse at large magnetic latitudes 
but are compressional near the equator. However, at the 
equator the transverse component of these waves increases 
relative to the compressional component as the degree of 
polarization decreases. Also, there is a substantial minority of 
events that are transverse at all magnetic latitudes, including 
the equator. A few of these latter events could be interpreted 
as ion cyclotron waves. In general, the waves tend to be 
strongly linear, have ellipticity magnitudes less than 0.3, and 
are skewed toward right-handedness, and the wave-nonnal 
angles are peaked near 90°. Polarization analysis indicates 
that a mixture of three wave modes is present, a left-handed 
(ICW) and a right-handed transverse mode and a compres- 
sional mode. For the transverse mode their wave magnetic 
field variations tend to be oriented perpendicular to the esti- 
mated field gradient, a pattern consistent with that of field- 
line resonances. If these transverse waves are field-line 
resonances, their standing wave pattern must be composed of 
at least two traveling waves, and therefore the interpretation 
of ellipticity and wave-normal angle from polarization anal- 
ysis is difficult [Denton et al., 1996], We note that for the 
transverse wave events plotted in Figure 9c, the events with 
negative ellipticity have larger power on average than those 
with positive ellipticity, which matches expectations if the 
waves are a superposition of two ICW modes. Despite the 
strong right-handed skew in the polarization analysis, there- 
fore, these waves may be composed mainly of compres- 
sional-mode waves and ICWs. Can the compressional and 
the transverse waves be related? The /3 H + (Table 2) in this 
wave environment can be very high, and Denton et al. [20 1 0] 
showed that in Earth’s plasma sheet boundary layer the 
compressional or transverse nature of the electromagnetic ion 
Bernstein wave mode can be strongly sensitive to /3 H +- 
Analysis of loss cone distributions in plasmas of highly var- 
iable but high /3 h + is needed to further delineate the role that 
these waves play in ion precipitation and in the direct gen- 
eration of field-line resonances at Mercury. 
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